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SUMMARY 


These  ranges  correspond  to  plasjnas  of  interest  for  magnetonydro- 
Only  monatomic  cesium  ions  are  considered  in  a 


fj  The  rate  of  three-body  recombination  between  cesium  seed  ions  and  elec¬ 
trons  in  an  argon  carrier  gas  is  calculated  for  free-electron  number  densities 
of  10^^  to  10^®  per  cubic  centimeter  and  electron  temperatures  of  5009^  to 
10,000°  K. 

dynamic  power  generation, 
collisional  approximation  where  radiative  transitions  are  ignored.  The  Byron 
method  is  applied  to  describe  collision-induced  transitions  between  atomic  en¬ 
ergy  levels  by  the  semiclassical  Gryzinski  cross  sections  and  the  principle  of 
detailed  balancing.  Recombination  is  limited  by  deexcitation  due  to  super¬ 
elastic  collisiona>-*^-  ^ 


Results  are  in  fair  agreement  with  those  of  a  more  approximate  study  of 
potassium  for  500°  to  2000°  K.  Discussion  of  the  results  as  related  to  conduc¬ 
tivity  in  magnetohydrodynamic  power  generation  and  comparison  with  recombina¬ 
tion  experiments  are  included.  The  computer  program  for  a  general  recombina¬ 
tion  study  of  an  atomic  ion  is  included  as  an  appendix.^ 


INTRODUCTION 

Recombination  Problem  in  the  Magnetohydrodynamic  Generator 

|The  requirements  of  magnetohydrodynamic  (MHD)  power  generation  set  lower 
limits  on  allowable  conductivity  values  for  gas  dynamic  systems  employing 
alkali-seeded  rare  gases.  The  electrical  conductivity  in  the  generator  duct  is 
a  function  of  the  free-electron  number  density  Ng,  which  must  be  maintained  by 
electron-atom  ionizing  collisions.  However,  if  the  electrons  initially  pro¬ 
duced  by  ionizing  collisions  recombine  with  positive  ions  in  a  time  that  i£| . G"  J 


*A  brief  summary  of  some  of  the  results  presented  in  this  report  was  given 
at  the  American  Physical  Society  Meeting,  Plasma  Physics  Section,  San  Diego, 
Calif.,  Nov.  2-6,  1963. 


f  short  compared  vith  a  characteristic  ionization  time^  volume  recombination  will 
severely  inhibit  attainment  of  the  conductivity  values  required  for  efficient 
generator  operation.)  A  recent  paper  predicts  high  recombination  rates  for  a 
potassium  ion  plasma  over  all  number  densities  in  an  electron  temperature  range 
from  500°  to  2000°  K  (ref.  l) .  TThe  object  of  the  present  report  is  to  estimate 
the  net  rate  of  electron-ion  recombination  in  a  plasma  of  argon  (ionization 
potential,  15.76  ev)  seeded  with  varying  percentages  of  cesium  vapor  (ioniza¬ 
tion  potential,  3.89  ev) .  These  gases  have  been  considered  a  promising 
carrier-seed  combination  for  magnetohydrodynamic  power  generation^) consequent¬ 
ly,  the  system  has  been  the  subject  of  extensive  study  (refs.  2  and  3). 


Selection  of  Electron  Capture  Mechanism 
F It  is  well  documented  that  three-body  capture 

Cs"^  +  e"  +  e"  ->  Cs*  +  e"  (K.E. ) 

is  the  initial  ^ep  of  the  neutralization  process  in  the  plasmas  of  Interest^ 
(refs.  1  to  4);lCs*  represents  an  excited  or  ground-state  seed  atomj  /This 
capture  step  is  a  radiationless  process  where  the  second  electron  carries  away 
the  excess  energy  as  kinetic  energy .J[  D'Angelo  performed  a  calculation  on  a-^  f  J' 
fully  ionized  gas  composed  only  of  protons  and  electrons  (ref.  4).  He  postu¬ 
lated  the  aforementioned  three-body  initial  capture  step  and  predicted  recom¬ 
bination  rates  in  agreement  with  several  sets  of  experimental  data,  but  his 
model  has  since  proven  inadequate  for  dense  plasmas  (Wg  >  lO^^  cm"^)  because 
he  ignored  deexcitation  collisions  (ref.  5).  Bates,  et  al.  (ref.  2)  and  the 
investigators  of  references  1  and  3  have  performed  more  general  studies  in 
which  they  have  included  this  effect. 

It  was  proposed  in  reference  6  that  the  neutral  carrier  atom  is  the  third 
body  for  certain  cesium  seeded  systems.  This  cannot  be  the  case  for  the  seeded 
argon  system  because  of  efficiency  requirements  for  energy  transfer  in  the 
three-body  encounter.  The  electron,  on  the  average,  can  transfer  only  the 
fraction  2mg/M  of  its  translational  energy  to  a  neutral  argon  atom  in  an 
elastic  collision  where  mg  and  M  are  the  electron  and  neutral  masses,  re¬ 
spectively.  Since  the  first  excited  state  of  argon  lies  11  electron  volts 
above  the  electronic  ground  state,  only  very  few  electrons  could  transfer  en¬ 
ergy  inelastlcally  to  the  atomic  third  body  for  the  6^  values  of  Interest. 

As  a  result,  only  electrons  in  the  extreme  low-energy  portion  of  the  Maxwellian 
distribution  can  successfully  engage  in  a  one-act  capture  process  where  an  atom 
serves  as  the  third  body.  On  the  other  hand,  electrons  can  exchange  any  amount 
of  kinetic  energy  in  their  collisions  with  one  another  and  thus  would  be  the 
most  effective  third  bodies  for  recombination. 


Byron  Approach 

Once  electrons  are  captured  into  excited  states  of  the  atom,  they  are  not 
considered  to  have  recombined  effectively  with  the  ion  since  they  may  have  a 
high  probability  of  being  reionized.  Byron,  et  al.  (refs.  1  and  3)  have  shown 
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that  the  rate-limiting  step  in  the  "chain"  of  three-body-recombination  pro¬ 
cesses  is  the  deexcitation  of  captured  electrons.  Their  model  utilizes  the 
rate-process  principle  that  the  slowest  step  in  a  rate  mechanism  is  the  limit¬ 
ing  step.  As  the  system  approaches  equilibrium,  the  rate  at  which  electronic 
levels  are  crossed  in  the  downward  direction  toward  the  ground  state  determines 
the  net  recombination  rate.  This  is  true  because  once  the  electron  has  reached 
the  ground  state  it  has  a  much  smaller  probability  of  being  reionized  than  it 
had  in  the  higher  states. 

The  authors  of  reference  3  have  shown  that  there  is  a  minimum  in  the  elec¬ 
tronic  deexcitation  rate  at  a  particular  atomic  energy  level.  The  location  of 
this  level  depends  strongly  on  the  temperature  of  the  free  electron  gas. 

In  their  application  of  the  rate-limiting  model  to  potassium  ion-electron 
recombination,  Byron,  et  al.  assumed  that  the  electronic  levels  above  and  below 
a  quantum  state  L*  are  continuous  bands  of  energy  levels  rather  than  discrete 
states.  They  did  not  list  the  states  included  in  their  study  nor  the  degener¬ 
acies  assigned  to  those  states.  In  addition,  they  approximated  the  slope  of 
the  Gryzinski  cross  section  by  a  linear  function  in  the  electron  temperature 
range ,  500°  to  2000°  K,  over  which  the  calculation  was  performed. 

The  present  report  uses  the  same  cross  section  expressions,  in  an  exact 
fashion,  with  the  rate-limiting  model  just  described,  and  applies  the  method  to 
cesium  recombination.  It  proposes  to  assess  the  accuracy  of  the  approximations 
given  in  reference  3  and  extends  the  range  of  the  recombination  coefficients  to 
temperatures  from  2000°  to  10,000°  K.  Any  quantum  effects  obscured  by  the 
continuum  approximation  of  Byron,  et  al.  should  appear  in  this  more  exact 
study.  Further,  the  information  contained  herein  is  designed  to  provide  a 
basis  for  calculating  recombination  in  high-pressure  plasmas  in  addition  to 
assessing  the  assumptions  of  the  Maxwellian  distribution  function  and  the  ne¬ 
glect  of  radiation. 

The  computer  program  for  treating  this  three-body  recombination  of  any 
monatomic  species  is  described  in  appendix  D  by  L.  U.  Albers.  This  program  is 
included  so  that  investigators  can  readily  perform  machine  calculations  on  sys¬ 
tems  of  interest. 


PRESENT  STUDY 


Scope 

friie  recombination  model  adopted  herein  is  applied  to  calculating  net  re¬ 
combination  rates  for  electron  number  densities  Ng  from  10^^  to  10^®  parti¬ 
cles  per  cubic  centimeter  and  electron  temperatures  .  0^  from  500°  to  10,000°  K 
for  a  conventional  range  of  cesium  seed  fractions  from  5x10“'^  to  10“^  (All 
symbols  are  defined  in  appendix  A. )  These  are  the  ranges  of  interest  for 
magnetohydrodynamic  generators.  Number  densities  below  10^^  per  cubic  centi¬ 
meter  will  not  provide  sufficient  conductivity  for  generator  operation,  while 
number  densities  of  ionized  seed  N+  (where  N+  =  Ng)  greater  than  lO^Q  per 
cubic  centimeter  correspond  to  neutral-carrier -gas  densities  of  at  least  10^0 


LI-  - 
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region,  Ng  greater  than  10^®  per  cubic  centimeter,  significant  overlap  of  cap¬ 
ture  orbitals  for  electrons  may  preclude  employment  of  the  adopted  recombina¬ 
tion  model,  since  collective  interactions  become  important.  With  regard  to 
electron  temperatures,  it  appears  that  the  most  promising  systems  may  make  use 
of  the  high  nonequilibrium  0g  values  from  2000°  to  8000°  K  (ref.  7)j  however, 
the  electrons  must  be  "Joule  heated"  from  the  static  gas  temperature,  and  this 
necessitates  a  recombination  study  at  6^  values  down  to  500°  K.  This  temper¬ 
ature  is  a  safe  lower  limit,  since  the  coldest  electrons  will  be  at  least  in 
themial  equilibrium  with  gas  atoms. 


Collisional  Approximation 

It  is  assumed  that  all  deexcitation  is  collisional  so  that  all  radiative 
deexcitations  can  be  neglected.  The  validity  of  this  assumption  is  discussed 
in  appendix  B,  which  indicates  that  radiation  from  excited  seed  atoms  is  negli¬ 
gible  down  to  Ng  values  of  about  10^'^  per  cubic  centimeter  for  all  electron 
temperatures.  Two-body  electron-ion  radiative  recombination  is  also  negligible 
in  the  Ng  range  considered.  Further,  it  is  assumed  that  the  sole  ionic  spe¬ 
cies  of  importance  is  the  cesium  atomic  ion;  that  is,  the  relative  abundances 
of  argon-cesium  ions  (ArCs"^)  and  argon  ions  are  very  small  compared  with  the 
cesium  ion  concentration.  This  is  a  good  approximation,  since  the  argon-cesium 
ion  (ArCs+)  is  bound  by  weak  polarization  energy  of  less  than  thermal  energy 
even  at  the  lowest  temperatures,  while  ionization  of  argon  will  have  a  rela¬ 
tively  long  relaxation  time  compared  with  cesium  ionization  (ref.  8).  Thus, 
argon  can  effectively  be  neglected;  preliminary  calculations  Indicate  that  this 
is  a  good  approximation  to  seed  fractions  of  2x10““^  for  temperatures  below 
8500°  K.  At  higher  electron  temperatures  there  may  be  some  ionization  of  ar¬ 
gon,  at  least  near  the  entrance  region  of  the  generator  duct  for  representative 
gas  residence  times  (ref.  9). 

Distribution  Function  Considerations 

Throughout  this  study,  it  is  assimied  that  a  Maxwellian  distribution  of 
free  electrons  is  preserved  in  the  plasma.  This  is  a  reasonable  assumption  for 
the  argon-cesium  seed  system,  as  can  be  shown  from  collision-frequency  consid¬ 
erations.  The  electron-electron  monoenergetic  collision  cross  section  for  mo¬ 
mentum  transfer  at  these  low  energies  (<1  ev)  is  relatively  high,  as  seen  even 
for  the  hard-sphere  low  estimate,  which  can  be  written  as  Qgg  =  SSS/w^  ang- 
strom^  (ref.  lO).  The  quantity  w  is  the  kinetic  energy  of  the  electron  in 
electron  volts.  The  related  collision  frequency  is  given  by  =  N^Q^gV^, 

where  Vg  is  the  relative  velocity  of  the  electrons.  The  frequency  of 
electron-neutral  encounters  is  given  by  Vg^^  =  ^n‘^n'''’en->  "wl^sre  Nj^  is  the 
neutral  argon  number  density,  Qg^^  is  the  elastic  electron-neutral  collision 
cross  section,  and  Vg^  is  the  electron-neutral  relative  velocity.  Even 
though  Ng  is  less  than  or  equal  to  10“^  N^  for  the  plasmas  of  interest,  the 
V  value  will  nevertheless  be  greater  than  because,  due  to  the  Ramsauer 

effect  (ref.  ll),  the  argon  atoms  are  virtually  transparent  to  electrons  in  the 
range  of  energy  from  0.1  to  1.0  electron  volt,  so  that  the  average  Qgj^  is 


only  3  angstroms^.  This  value  is  the  average  value  of  from  0.1  to 

1.0  electron  volt  for  Maxwellian  distributions  of  electrons  at  each  electron 
temperature.  Elastic  collisions  between  electrons  and  seed  atoms  are  assumed 
to  have  maximum  cross  sections  of  the  order  of  the  value.  This  latter 

assumption  is  well  founded  since  the  (pessimistic)  minimum  (at  1  ev)  is 

656  angstroms^^  while  the  maximum  calculated  for  cesium  neutrals  is 

600  angstroms^  (averaged  over  the  Maxwellian  distribution  at  0.5  ev^  ref.  12). 
Consequently^  the  electron  collision  frequency  is  certainly  greater  than  either 
electron-neutral  (carrier  or  seed)  collision  frequency.  The  assumption  of  a 
Maxwellian  electron  distribution  should  therefore  be  valid. 


DETAILS  OE  CALCULATION 


The  deexcitation  rate  R.  ,  is  defined  herein  as  the  rate  at  which  elec- 

dex 

trons  in  discrete  energy  states  K  above  a  state  L"^  can  pass  downward  either 
to  L*  or  past  that  state  to  all  states  L  below  it^  across  the  energy  gap 
L*,  where  K*  =  L*  +  1.  It  is  assimied  that  all  the  electronic  states  are 
in  equilibrium  with  the  free  electrons  at  the  electron  temperature.  This  is 
not  quite  true  for  the  lower  states^  whose  populations  are  less  than  the  equi¬ 
librium  values  because  they  require  a  longer  time  to  attain  equilibrium  popula¬ 
tions.  To  account  for  these  and  other  nonequilibrium  effects^  a  temperature- 
dependent  factor  X  introduced  to  relate  the  minimum  deexcitation  rate  to 
the  recombination  rate. 


For  a  hypothetical  nonequilibrium  situation  (l)  with  no  reexcitation  oc¬ 
curring  during  the  period  of  deexcitation  and  (2)  state  populations  above  the 
limiting  gap  equal  to  equilibrium  values^  y  would  be  unity^  and  the  net  recom¬ 
bination  rate  would  be  equal  to  the  deexcitation  rate  for  the  slowest  or  rate¬ 
determining  transition.  Clearly^  at  equilibrium  where  =  ^ex 

gaps^  there  is  no  net  recombination  and  y  ^must  be  zero-  The  calculation  of 
recombination  rate^  therefore^  consists  of  determining  the  minimum  deexcitation 
rate  and  the  value  of  y  for  the  range  of  electron  temperatures  and  densities 
of  interest. 


Collision  Cross  Sections  and  Deexcitation 

For  a  one-electron  outer  shelly  such  as  that  of  the  singly  excited  or 
ground  state  cesium  atom^  the  ionization  potential  (or  binding  energy)  is  equal 
to  the  kinetic  energy  of  the  bound  electron.  Gryzinski  (ref.  13)  has  derived 
cross  section  expressions  for  the  inelastic  collision  of  an  electron  with  an 
atom  by  treating-  such  a  collision  as  a  classical  encounter  between  the  free 
electron  and  the  bound  electron.  Although  these  expressions  predict  a  low 
value  and  the  wrong  energy  dependence  of  the  ionization  cross  section  at  the 
threshold  (aEl/2  instead  of  the  experimentally  observed  oc  E  (ref.  14))^  they 
should  be  satisfactory  for  calculating  the  recombination  rate  within  a  factor 
of  two  to  fo-ur  (ref.  15).  Experimental  error  in  recombination-rate  measure¬ 
ments  is  somewhat  less  (refs.  5  and  16)^  so  that  discrepancies  between  pre¬ 
dicted  and  observed  rates  can  be  explained^  at  leasts  in  part^  as  due  to  in¬ 
adequate  cross  sections. 
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The  main  deexcitation  route  in  the  range  of  interest  for  Ng  and  6^  is 
through  superelastic  collisions  whereby  the  incident  (free)  electron  receives 
the  energy  of  atomic  deexcitation  from  the  bound  electron  as  increased  kinetic 
energy  (refs.  1,  2,  and  16). 

In  the  Byron  approach^  the  deexcitation  rate  is  not  computed  directly 
through  the  deexcitation  cross  section  but  must  be  obtained  from  the  excitation 
rate  by  using  the  statistical  mechanical  principle  of  detailed  balancing 
(ref.  17).  The  inelastic  excitation  cross  section  necessary  to  calculate  the 
excitation  rate  is  discussed  in  the  following  section.  The  Gryzinski  expres¬ 
sions  have  been  remarkably  successful  (in  view  of  their  derivation  on  a  purely 
semiclasslcal  basis)  in  predicting  collisional  excitation  probabilities  from 
ground  to  lower  excited  states  for  some  atoms. 


Energy  Levels  of  Cesium  Atom  Quantum  Degeneracy 

The  first  70  observed  spectral  levels  of  the  cesium  atom  selected  for  the 
calculation  are  taken  from  reference  18.  Because  of  pressure  and  Doppler 
broadening  (ref.  l)  and  external  magnetic -fie Id  effects  (B  may  be  as  high  as 
20,000  gauss),  it  was  thought  reasonable  to  assume  that  energy  states  separated 
by  30  wave  numbers  or  less  are  effectively  degenerate  in  energy.  The  47  elec¬ 
tronic  "uncollapsed"  states  are  listed  in  table  I  with  their  binding  energies 
E-,  .  The  n  and  I  values  are  the  principal  and  secondary  (orbital)  quantum 
numbers,  whereas  Sp  =  l/2  is  the  spin  of  one  unpaired  electron  in  units 
of  ft. 


The  total  degeneracy  of  an  electronic  energy  level  o)  is  given.  Spin- 
orbit  coupling  removes  the  spin  degeneracy,  giving  rise  to  the  spin-orbit  quan¬ 
tum  number  J  =  Z  ±  Sp. 

Eighteen  of  the  listed  levels  are  nondegenerate  in  J  value.  The  remain¬ 
ing  29  levels  denoted  by  asterisks  are  degenerate  levels.  Seven  of  these 
states  degenerate  in  J  are  s  states  (Z  =  O)  where  the  spin  degeneracy  can 
only  be  removed  in  an  external  magnetic  field.  The  other  22  states  are  p,  d, 
f,  g,  and  h  states  (Z  =  1  to  5)  where  Z-s  coupling  is  not  strong  enough  to 
split  the  levels  by  more  than  30  wave  numbers  (cm“^)  in  energy. 


Excitation  Cross  Section 


The  cross  section  (cm^)  for  inelastic  collision  "between  an  incident  elec¬ 
tron  with  energy  E^  and  an  atomic  electron  with  kinetic  energy  (or  binding 
energy)  E^^  occurring  with  an  energy  loss  greater  than  AE  is  (ref.  13) 


(la) 


for  Z!^  +  Ej  <  Eg 
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and 


In  both  expressions 


for  AE  +  >  Eg 


(4.8x10~^Q)^ 

(1.6X10"^^)^ 


=  6.53x10“!^  (cm2)(ev2) 


(lb) 


where  e  is  the  electronic  charge. 

The  mean  product  of  cross  section  and  velocity  (cm^/sec),  integrated  over 
the  Maxwellian  distribution^  for  the  excitation  of  an  electron  from  an  energy 
level  L  up  to  an  energy  level  K  is  (ref.  19^  pp.  677-699) 


(''Wl-k))  =-g- 


[%,,K 


f 


G(E2,Ei,ABl,k) 


G(E2;E2_>AEj^^K+i) 


,K+1 


F(E;L>^2)^2  <^2  (2) 


where  G(E2^E]_,AE)  is  the  fimction  in  braces  of  equation  (la)  or  (ib)^ 

F(E^,E2)  =  [e2/(E^  +  E2)]^/^  g“^2A®e^  a.nd  Cq  is  the  normalization  factor 
for  the  Maxwellian  distribution  f-unction^  (2jtk0g)^/^(mg)^/^. 
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An  electron  that  is  excited  above  the  level  K,  (AE  >  ^  but  not  as 
high  as  the  level  K  +  1,  (AE  <  AE^  k+i)  •»  reach  the  level  K  by  emitting 
radiation.  This  is  a  physical  basis  for  the  tvo  lower  limits  in  the  integral 
expression  for  the  mean  excitation  cross  section. 

The  appropriate  should  be  employed  for  the  two  sets  of  conditions 
mentioned  in  equations  (la)  and  (lb). 


Excitation  Rate  and  Detailed  Balancing 

The  corresponding  rate  at  which  electrons  can  be  excited  from  the  state  L 
to  state  K  is 

The  rates  of  the  inverse  processes^  excitation  and  de excitation^  are  equal 
provided  that  the  sum  of  the  frequencies  for  super elastic  collisions  from  each 
state  K  to  all  possible  states  L  is  larger  than  the  sum  of  the  Einstein  co¬ 
efficients^  which  are  the  transition  probabilities  for  radiative  emission  from 
the  upper  state.  This  condition  defines  detailed  balancing^  according  to 
which ^  if  a  Maxwellian  distribution  of  electrons  is  assumed^  the  deexcitation 
rate  ^(3_ex  between  any  two  states  is  equal  to  the  excitation  rate  ^ex* 
deexcitation  rate  from  state  K  to  state  L  can  be  written 


'dex 


AEp^^L/ke, 


((Ql^k)t) 


(4) 


Continuum 


K  =  K*  +  1 


K  =  K*^ 

»  Sample  gap 
L  =  L"  =  4 

L  =  3 
L  =  2 


Ground  _ _ .  _ , 

state 

Figure  1.  -  Schematic  energy  level  diagram  for  cesium  atom  that  shows 
sample  gap. 


where 


_  K\  ^,LA®e 


is 


the  ratio  of  the  deexcltatlon  colli¬ 
sion  coefficient  to  the  excitation 
collision  coefficient.  In  equa¬ 
tion  (4)  a>]^  and  are  the  total 

degeneracies  of  the  levels  L  and 

K  =  +  !)• 

J 


The  rate  of  deexcltatlon  R^ex 
across  a  gap  K*  L*  is  the  svim  of 
the  rates  at  which  K*  and  all  states 
K  above  the  state  can  be  de- 

excited  by  superelastic  collision 
either  past  the  state  L*  to  all 
states  L  below  L*  or  to  the  state 
L*.  Figure  1  schematically  displays 
the  arrangement  of  energy  levels  and 
gaps  for  cesium  below  the  electron 
continuum^  which  is  indicated  by  the 
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solid  dark  area.  The  rate  can  be  written 

K=N  L=L* 

%ex  =  ^  ^  ^dex 

K=K*  L=1 


(5) 


where  all  states  K  are  above  K*  and  all  states  L  are  below  L*  as  indi¬ 
cated  in  the  figure,  and  N,  the  total  number  of  states,  is  equal  to  47. 


In  all  cases,  the  population  (cm"^)  of  the  state  K  is  given  by  the 
Saha  equation  since,  as  mentioned,  the  bound  states  are  assumed  to  be  in  equi¬ 
librium  with  the  free  electrons.  Each  population  is  given  by 


'K 


where  is  the  binding  energy  of  level  K  (its  previously  introduced  E3_ 

value),  cDg  is  2  for  the  free  electron,  and  cd^  is  1  for  the  degeneracy  of  the 
cesium  ion. 


Estimation  of  y.  -  The  value  of  y  in  the  relation  =  r(R.q  )  . 

IS  an  indication  of  the  magnitude  of  nonequilibrium  effects  (refs.  1  and  3) 
and  is  a  function  of  electron  temperature.  This  factor  must  account  for  the 
fraction  of  the  minimum  deexcitation  rate  that  is  balanced  by  excitation  as  the 
system  approaches  equilibrium,  as  well  as  for  nonequilibrium  populations. 

The  deexcitation  fraction  as  well  as  the  excitation  and  deexcitation 
cross  sections  for  important  transitions  across  the  gap  are  calculated  by  com¬ 
paring  relaxation  times  required  to  attain  equilibrium  for  states  above  and 
below  the  minimum.  The  ratio  of  these  cross  sections  contains  the  factor 
,  .  ,  AEr  v-/ke„ 

(a^/a^)e  ’  where  is  the  energy  difference  and  the  co  values 

are  the  degeneracies  for  the  significant  transitions. 

At  the  lower  electron  temperatures,  where  the  rate -controlling  gap  is  lo¬ 
cated  between  high  quantum  states,  neither  capture  rates  nor  state  populations 
differ  much  for  adjoining  states  above  and  below  the  gap,  consequently,  re¬ 
excitation  is  Important.  The  energy  differences  for  important  transitions  are 
less  than  kOg,  so  that  the  excitation  cross  section  is  usually  comparable  to 
the  deexcitation  cross  section.  Therefore,  for  low  electron  temperatures, 

500°  to  3000°  K,  the  deexcitation  fraction  must  first  be  determined  by  detailed 
examination  of  the  significant  excitation  cross  sections,  and  then  the  popula¬ 
tions  of  the  states  above  the  gap  are  corrected  for  their  departure  from  equi¬ 
librium  values  to  obtain  the  y  factor. 

At  high  temperatures,  where  the  (6p  6s)  gap  controls  the  rate  process, 

capture  rates  are  very  different  for  states  above  the  gap  and  the  ground  state, 
so  that  the  deexcitation  fraction  is  near  unity  and  y  must  account  only  for 
the  Inaccurate  populations  of  the  6p  state  and  several  states  above  it. 
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In  either  case  any  inaccuracy  involved  in  estimating  y  is  veil  vithin 
the  error  due  to  the  collision  cross  sections.  These  considerations  yield  es¬ 
timates  for  r  of  1/3  through  2000°  K,  l/2  from  3000°  to  5000°  K,  and  2/3 
for  5000°  to  10,000°  K. 


RESULTS  AHD  DISCUSSION 


E*  = 
ure 


The  plot  of  sec"!)  against  the  mean  gap  energy, 

Ej^  +  ^  (Ej^  -  Ejr)  calculated  from  equations  (4)  and  (5),  is  given  in  fig- 


2  for 


=  1000" 


K  and  N, 


=  10^3  per  cubic  centimeter.  The  values  for 

these  energy  gaps  are  given  in 
table  II.  The  dashed  curve  is 
a  similar  plot  from  reference  1 
for  potassium  and  shovs  a  more 
regular  behavior  than  the 
points  obtained  in  the  present 
study.  It  should  be  noted  that 
Byron,  et  al.  (ref.  l)  assumed 
that  the  levels  above,  as  veil 
as  belov,  L*  are  continuous 
rather  than  discrete  and  also 
that  the  slope  of  the  Gryz inski 
cross  section  can  be  approxi¬ 
mated  by  a  linear  function  in 


Figure  2.  -  Deexcitation  rate  as  function  of  mean  gap  energy.  Eiectron 
number  density,  10^^  per  cubic  centimeter;  free-eiectron  temperature, 
1000“  K. 


the  range  of  6g  from  500°  to 
2000°  K.  There  is  more  quantum 
fine  structture  in  the  results 
of  this  study  although  the  min¬ 
imum  R(3^gjj(vhich  gives  R^-ec^ 
is  not  too  different  in  E* 
value  from  reference  1.  This 
agreement  is  expected  since 
these  tvo  alkalis  have  qualita¬ 
tively  similar  arrangements  of 
atomic  energy  levels. 

The  locus  of  minimums  in 

figure  2  denotes  the  gaps 

K*  ->  L*  that  have  the  lovest 

R.  values  in  the  various 
hex  „ 
ranges  of  ^ 


=  1000°  K. 


E*  at 

The  recombination  rate  is  one- 
third  the  value  of  the  minimum 


of  the  locus  curve,  vhich  occurs  in  this  case  for  the  (l2s,  8f  lOd)  gap  at 


E*  =  0.23  electron  volt. 


FigfH’e  3  shovs  the  recombination  rate  Rpec 
temperature  0g  in  the  range  500°  to  10,000°  K. 
of  the  rate-limiting  mean  gap  energy 


as  a  function  of  electron 
It  is  clear  that  the  location 
E*  is  depressed  vith  increasing  tempera- 
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Free-electron  temperature,  eg,  °K 


Figure  3.  -  Three-body  coliisional  recombination  coefficient  as  function  of 
free-electron  temperature  for  range  of  electron  number  density  10^^  to  10^® 
per  cubic  centimrter. 


ture.  It  should  he  noted 
that  the  (6p  6s)  energy 

difference  results  in  a 
relatively  long  time  for 
the  ground  state  to  come 
to  equilibrium  so  the  de¬ 
excitation  is  only 
slightly  balanced  by  ex¬ 
citation.  At  tempera¬ 
tures  above  10^000^ 
two-body  radiative  recom¬ 
bination  becomes  impor¬ 
tant  along  with  quantum 
effects^  so  that  the 
purely  coliisional  form 

=  \ec/^e 
longer  adequate. 

The  temperature  de¬ 
pendence  in  the  range 
500°  to  2000°  K  is  much 
the  same  as  the  dashed 
curve  shown  for  the  work 
of  Byron^  et  al.  (ref.  l) 
in  figure  3.  This  tem¬ 
perature  coefficient  a* 
is  similar  to  the  a  of 
reference  2^  which  is  a 
treatment  of  a  quite  dif¬ 
ferent  nature  from  the 
present  calculation.  Ad¬ 
ditional  calculations 


that  were  performed  with 

59  discrete  energy  levels  showed  that  *'collapse"  of  energy  levels  separated  by 
small  gaps^  which  produces  new  cd  values^  changes  the  calculated  rates  by  only 
several  percent  at  most.  The  temperature  dependence  of  a*  goes  from  0“^  in 
the  range  500^  to  3000°  K  to  0“^/^  in  the  range  of  0^  from  4000°  to 

10^000°  K  compared  with  the  0^^/^  dependence  of  reference  5. 


The  results  shown  in  figure  3  are  in  good  agreement  (within  a  factor  of  2) 
with  the  experimental  findings  of  reference  20  for  a  cesium  plasma  with 
values  from  IxlO^^  to  6xl0l^  per  cubic  centimeter  and  0^  values  from  0.11  to 
0.26  electron  volt.  They  also  agree  closely  with  the  data  of  reference  21 
taken  for  a  pure  thermal  cesium  plasma  in  the  ranges  of  of  5x10^^  to 

10l3  cubic  centimeter  and  0^  of  2000°  to  2500°  K.  Since  both  radiative 
deexcitation  and  trapping'^of  resonance  radiation  occur  in  these  plasmas  for 
values  of  interest^  it  is  pointless  to  offer  a  refined  interpretation  of  the 
difference  between  experiment  and  theory. 


The  theoretical  results  are  also  in  good  agreement  with  the  low  electron 
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temperature  experimental  data  of  reference  22  for  an  argon  plasma  seeded  with 
cesium,  which  is  in  excellent  agreement  with  the  theoretical  value  of  refer- 
ence  5. 


The  recombination  results  for  a  modified  D'Angelo  approach  (ref.  4)  are 
shown  for  comparison  in  figure  4.  This  supplementary  calculation  was  done  in 


a  manner  quite  different  from  the  .Byron  approach  and  is  outlined  in  appendix  C. 


COWCLUDIWG  KEMAEKS 

|"For  conventional  output  power  densities  of  the  magnetohydrodynamic  gener¬ 
ator,  recombination  could  be  a  problem,  especially  below  free-electron  tempera¬ 
tures  of  3000°  K  because  of  the  strong  inverse  temperature  dependence  of  the^ 
net  recombination  rate.  However,  as  others  have  suggested -(Byron, .et  al. )  with 
the  use  of  electron  beams  to  "kick  up"  the  electron  temperature  initially,  the 
induced  field  in  the  generator  may  then  provide  "sustained"  electron  heating  to 
keep  the  electrons  at  an  elevated  temperature,  thus  avoiding  the  harmful  recom¬ 
bination  losses.  With  the  "hot"  electrons  (electron  temperature  »  gas  temper- 
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ature)  the  cesiian  seeded  argon  system  remains  the  most  promising  one  for  ob¬ 
taining  high  concjg^ivity  at  optimum  operating  conditions  in  view  of  the  fact 
that  the  net  recombination  rates  calculated  for  plasmas  of  interest  are  not 
prohibitive .  j  ^ 

Although  the  acciuracy  of  the  method  is  limited  by  the  cross  section  ex¬ 
pressions  that  underestimate  the  excitation  rate  and,  consequently,  the  de¬ 
excitation  rate,  it  may  be  possible  to  ascribe  differences  between  experiment 
and  theory  to  radiation  effects  in  the  plasma.  For  example,  if  the  observed 
recombination  rate  is  significantly  less  than  the  theoretical  rate ,  the  differ¬ 
ence  may  be  attributed  to  radiation  trapping.  The  method  does  provide  a 
straightforward  approach  for  estimating  net  recombination  and  has  been  remark¬ 
ably  successful  thus  fax. 

The  computer  program  for  the  recombination  calculation  is  included  in  ap¬ 
pendix  D  with  explanatory  remarks.  For  any  number  densities  or  electron  tem- 
perat-ures,  the  only  additional  data  required  are  the  energy  levels,  values, 
and  corresponding  degeneracy  values  co  for  the  atom  of  interest.  The  y  fac¬ 
tor  can  be  estimated  as  described  in  the  text. 


Lewis  Research  Center 

Rational  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  July  27,  1964 
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APPEKDIX  A 


SYMBOLS 

Einstein  coefficient  for  emission  Between  states  K  and  sec“^ 
normalization  factor  for  Maxwellian  distribution 
inelastic  energy  loss  for  collision^  ev 
binding  energy  of  level  ^lK-> 

kinetic  energy  of  bound  electron  (ionization  potential  of  atomic  state 
ev 

kinetic  energy  of  incident  electron^  ev 
mean  binding  energy  of  quantum  gap^  ev 
unit  electronic  charge^  4.8x10”^^  esu 
Planck ^s  constant^  4.1x10"^^  (ev)(sec) 
h/2jt 

spin  orbit  quantum  number^  Z  ±  Sp 

electronic  energy  level  above  gap;  K*  the  electronic  energy  level 
immediately  above  gap 

Boltzmann^s  constant^  8.6x10“^  (ev)(^K"^) 

electronic  energy  level  below  gap;  L*  the  electronic  energy  level 
immediately  below  gap 

secondary  (orbital  angular  momentum)  quantum  number 

neutral  mass^  6.8x10“^^  g  (argon) ^  2.2x10“^^  g  (cesium) 

electronic  mass^  9.1x10"^^  g 

free-electron  nimiber  density^  cm“^ 

population  of  state  cm 

population  of  state  L,  cm"^ 

'Z 

ion  number  density,  cm 


principal  quantum  number 


’ee 


'en 


R 


dex 


R 


'ex 


rec 


ee 


en 


a ' 


ee 


en 


electron-electron  collision  cross  section  for  momentum  transfer,  cm^ 
electron-neutral  collision  cross  section  for  momentum  transfer,  cm^ 
deexcitation  rate  of  bound  electrons,  (cm“^) (sec"^) 
excitation  rate  of  bound  electrons,  (cm"^) (sec“^) 
net  recombination  rate,  (cm"^)(sec“l) 

total  spin  q^uantum  nimiber,  l/2  for  all  cesium  atom  states 
electron-electron  relative  velocity,  (cm)(sec"^) 
electron-neutral  relative  velocity,  (cm)(sec“^) 
free -electron  kinetic  energy,  ev 

collisional  recombination  coefficient,  (cm®)(sec“l) 

adjustment  factor  to  account  for  absence  of  complete  bound-free  equilib¬ 
rium  and  reexcitation  effects 

temperat-ure  of  free  electrons,  °K 

cross  section  constant,  6.53x10"^'^  (cm^)(ev^) 

electron-electron  collision  frequency  for  momentum  transfer,  sec"l 
electron-neutral  collision  frequency  for  momentum  transfer,  sec“^ 
statistical  weight  (total  degeneracy)  of  an  electronic  level  K, 

E(2J  +  1) 

J 
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APPEKDIX  B 


IffiGLECT  OF  RADIATION 

The  extent  to  which  radiative  emission  can  be  neglected  as  a  deexcitation 
mechanism  can  be  assessed  from  table  III^  which  lists  the  superelastic  colli¬ 
sion  frequencies  (sec”l)  for  all  deexcitation  collisions  from  a 

state  K.  Each  frequency  is  a  reciprocal  of  the  superelastic  relaxation  time. 
This  quantity  is  the  total  frequency  for  deexcitation  across  a  gap 

K*  L*  from  a  state  K  for  Ng  of  10^^  per  cubic  centimeter  and  6q  of 
5000°  K.  In  comparison^  the  oscillator  strengths  for  cesium  radiation  emission 
are  relatively  small  for  transitions  between  excited  states^  while  transitions 
to  the  ground  state  have  a  maximum  frequency  Aj^^  of  3x10''  per  second 
(ref.  23).  Most  of  the  Einstein  coefficients  Aj^  are  at  least  an  order  of 
magnitude  less  than  this  maximum.  Since  Xv^gj^  is  much  greater  than  XAjq^ 

it  appears  that  radiative  deexcitation  will  only  become  significant  for  Rg 
less  than  3x10^'^  per  cubic  centimeter  at  even  the  lowest  electron  temperatures. 
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APPEMDIX  C 


MODIFIED  D'ANGELO  MODEL 

A  separate  approach  to  the  recomhination  rate  was  made  in  this  study  with 
a  modification  of  the  D'Angelo  model.  Reference  4  outlines  a  treatment  whereby 
electrons  are  traced  from  the  continuum  after  being  captured  in  a  three-body 
process  at  a  rate  I^capJ  however,  the  only  downward  process  Included  between 
electronic  levels  was  radiative  transition  and  the  only  upward  process  included 
was  collisional  ionization.  This  study  includes  superelastic  collisions  and 
excitation  collisions.  The  inclusion  of  these  processes  leads  to  a  divergent 
recombination  rate  if  Bohr -Thomson  cross  sections  (see  ref.  16)  are  employed. 
This  occurs  because  capture  is  the  limiting  step  with  the  inclusion  of  colli- 
sional  processes  in  the  modified  model.  The  capture  rate  naturally  diverges 
for  increasing  orbital  size  and  decreasing  ionization  potential  or  binding 
energy;  however,  if  the  bound  state  distribution  is  fitted  to  the  Maxwellian 
distribution  in  the  manner  of  a  continuous  function,  as  suggested  by  John  E. 
Heighway  of  Lewis,  a  cutoff  is  derived  for  the  botind  states  at  3/2  (k0  ),  the 
average  of  the  free-electron  energies.  The  recombination  rate  is  then  calcu¬ 
lated  by  summing  over  all  states  below  the  cutoff. 


17 


APPENDIX  D 


COMPUTER  PROGRAM  FOR  STUDYING  RECOMBINATION  OF  AN  ATCMIC  ION 

by  Lynn  U.  Albers 

The  computer  program  was  ‘wri't’ten  in  FORTRAN  IV  compiler  language  for  exe¬ 
cution  on  an  IBM  7094  computer  under  control  of  the  Lewis  IBSYS  Monitor  System. 
It  consists  of  a  main  program,  an  integration  subroutine  VQEXI,  and  a  function 
subprogram  F. 

The  purpose  of  the  main  program  is  to  compute  the  following  three  sums  for 
all  (L*,K*)  gaps  where  K*  =  L*  +  1  and  L*  varies  from  1  to  the  number  of 
states  N: 

L*  N+1 

RDEX  =  ZZ.  Fi(L,K)  (Dl) 

L=1  K=K* 

L*  N+1 

F2SUM  =  ZZ.  F2(L,K)  (D2) 

L=1  K=K* 

L* 

DEXF  =  ^  F2(L,K*)  (D3) 

L=1 


where 

Fi(L,K)  =  Nj^F2(L,K)  (D4) 

and 

F^iL,K)  =  -  Ek)  -  -  %+i)]  if  K  <  B 

=  P3I(El,El)  if  K  =  N  +  1  (D5) 

are  eq.uations  related  to  equations  (2)  and  (4)  in  the  text.  The  parameter  P3 
is  one  of  three  parameters  defined  by 


(D7) 
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by  using  the  notation  of  equations  (la),  (lb),  and  (2).  In  this  subroutine, 
the  whole  integrand  is  denoted  by  F(E2)  with  and  Z®  being  available  in 

COMMON.  Since  some  factors  in  the  integrand  are  likely  to  be  smaller  than  the 
10-38  limit  for  single-precision  floating-point  numbers,  precaution  is  taken  to 
compute  a  BIAS  for  the  exponent,  which  is  stored  in  COMMON  and  used  to  bias 
the  argument  of  the  exponential  function  every  time  the  integrand  is  computed 
by  means  of  the  subprogram  F.  A  compensating  scale  factor  SCAL  is  included  in 
the  P3  formula  for  PARS.  The  integration  is  performed  by  Simpson's  Rule  in 
the  routine  VQEXI.  The  integrand  has  the  character  of  varying  widely  in  a  very 
narrow  domain  at  the  left  of  the  interval  and  being  negligible  elsewhere.  The 
VQEXI  routine  searches  for  this  narrow  interval  where  the  integrand  is  greater 
than  10"5  times  F(AE)  and  uses  32,  64,  or  128  integration  steps  to  evaluate 
the  integral. 

The  Grysinski  expressions  for  AE  =  Eg  give  falsely  negative  integrands, 
which  might  throw  some  doubt  on  end-point  formulas.  However,  even  internal 
formulas  like  Gaussian  quadrature  formulas  gave  the  same  results,  so  no  error 
is  involved  in  the  method  used.  The  logic  of  the  VQEXI  routine  is  simple 
enough  that  it  could  be  replaced  by  some  other  process  if  further  improvement 
in  execution  time  is  desired. 

After  the  reading  in  of  energy  levels  and  their  identification  numbers,  as 
well  as  conditional  input  of  modified  omegas,  all  parameters  independent  of 
temperature  and  electron  number  density  are  computed.  The  remainder  of  the 
program  accomplishes  an  iterative  input  of  a  temperatiire  and  an  electron  den¬ 
sity  and  a  calculation  of  an  output  of  the  s-ums  defined  in  equations  (Dl), 

(D2),  and  (D3). 

An  earlier  version  of  the  program  computed  the  sxmi  for  each  gap  separate¬ 
ly.  Since  for  gaps  near  the  middle  this  would  involve  the  sum  of  many  terms, 
apd  the  computations  for  adjacent  gaps  overlap  considerably,  the  program  was 
modified  to  compute  the  sum  for  one  gap  from  the  simi  for  the  preceding  gap  by 
making  proper  additions  and  subtractions.  When  this  was  first  done,  each 
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addend  was  computed  as  needed,  once  when  added,  and  again  when  subtracted.  By- 
storing  each  addend  the  first  time  it  was  computed,  the  execution  time  could  be 
halved.  The  disadvantage  of  proceeding  by  augmentation  and  deletion  is  that 
sums  near  the  end  of  the  process,  though  composed  of  a  relatively  small  number 
of  terms,  are  masked  by  the  round-off  error  attendant  on  adding  and  subtracting 
large  terms.  Hence,  a  provision  is  made  for  recomputing  gap  sums  from  the  low 
end  conditionally  for  a  number  of  steps. 

If  one  wishes  to  shorten  execution  time  by  neglecting  terms  between  widely 
separated  states,  one  can  by  changing  LKB,  limit  LL  to  being  greater  than 
L*  -  KLD,  and  limit  KK  to  being  less  than  +  LKD.  This,  in  effect,  trun¬ 
cates  the  sums  in  equations  (Dl)  and  (D2). 

V 

As  RDEX  and  F2SUM  are  computed  for  each  gap,  starting  at  the  low  end,  they 
are  listed  alongside  of  the  labels  for  L*  and  K*,  the  L*  and  K*  indices, 
and  the  population  of  the  upper  stat'e.  Then  if  HBW  is  positive,  gap  sums 

are  recomputed  and  listed  starting  from  the  high  end  for  the  required  number  of 
steps  by  using  the  F2  terms  stored  earlier.  The  table  TAB  has  five  rows  of 
99  elements  each.  F2SUM  and  RDEX  in  the  forward  sweep  are  stored  in  rows  1 
and  2.  F2SUM  and  RDEX  recomputed  in  the  backward  sweep  are  stored  in  rows  3 
and  4.  Finally,  the  DEXF  is  computed  and  stored  in  row  5.  If  a  TABLE  1  value 
is  negative,  this  indicates  round-off  error  for  this  and  all  higher  gaps,  and 
TABLES  1  and  2  are  over-written  by  the  more  acciurate  TABLES  3  and  4  from  this 
gap  to  the  top  gap.  It  may  be  necessary  to  recompute  the  tables  for  a  large 
number  of  states  from  the  high  end  at  lower  electron  temperatures  (e.g., 

WBW  =  43  at  ©g  =  1000°  K)  while  at  higher  temperatures,  that  is,  5000°  K  and 
above,  no  recomputation  is  necesseiry. 

Once  this  table  has  been  assembled,  it  is  read  out  in  two  sections  called 
table  II  and  table  III.  Table  II  lists  F2SUM,  RDEX,  and  the  average  ionization 
potential  for  the  gap  in  sets  of  twenty  in  such  a  -way  that  each  sum  is  brack¬ 
eted  by  the  labels  and  of  the  bounding  states  of  its  gap.  Table  III 

lists  DEXF  and  RDEX  in  corresponding  sets  of  twenty.  After  this  is  done,  a 
new  (0g,Wg)  pair  is  read  in  until  all  cases  are  treated. 

Following  is  an  outline  of  the  main  program  using  FORTRAN  symbols.  A 
dictionary  of  these  symbols  with  their  meanings  is  given  after  the  outline,  and 
finally  a  listing  of  (l)  MAIN,  (2)  VQEXI,  (3)  F,  and  (4)  a  typical  set  of  data. 

OUTLINE  OF  MAIN 

I.  Specify  constants  PI,  K,  H,  SIGZ,  and  ME. 

II.  Specify  atomic  data  GSIP,  NSTAT,  WE,  and  LKD. 

III.  Read  in  arrays  LAB,  EN,  and  OMEGA, 
rv.  List  OMEGA  array,  GSIP,  WE,  and  EN  array. 

V.  Provide  for  states  beyond  N^^. 
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VI.  Compute  IP  array  (loop  l)  and  list  it. 

VII.  Read  in  KE,  THETA,  HB¥  and  list  them. 

VIII.  Calculate  KT,  KTHET,  PARI,  PAR2,  and  IHF. 

IX.  Compute  M  array  (loop  2)  and  list  it. 

X.  Calculate  BIAS  and  SCAB. 

XI.  Skip  to  next  page  and  print  headings  at  top.  Set  F2SUM  and  RDEX  to  zero. 

XII.  Make  forward  sweep  calculating  and  storing  F2SUM  and  RDEX  for  each  gap  in 
TAB  matrix,  and  listing  them  and  associated  data  line  hy  line  (loop  3). 
LSTAR  runs  from  1  to  NSTAT. 

A.  Compute  KSTAR,  KMAX,  IMIN,  and  El. 

B.  Augment  siims  hy  proper  terms  for  this  gap  (loop  4)KK  =  KSTAR  to  KMAX. 

1.  Compute  DELE,  PARS  for  (L*,KK)  pair. 

2.  Do  first  integral  (see  eq..  (D5))  and  store  in  TERM. 

3.  Conditionally  compute  the  second  integral  for  (L*,KK  +  l)  pair 
and  subtract  it  from  TERM  to  get  F2  value. 

4.  Increment  F2SUM,  RDEX.  Store  F2  element  in  (L*,KK)  slot. 

C.  Subtract  proper  terms  from  sums  for  this  gap  (loop  5)  LL  =  LMIH  to 
L*  -  1. 

1.  Obtain  F2(LL,L*)  from  storage. 

2.  Decrement  F2SUM  and  RDEX. 

D.  Store  final  F2SUM  and  RDEX  in  Tables  1  and  2  in  L*  slot. 

E.  List  LAb(L*),  LAB(K:^),  L*,  K*,  F2SUM,  RDEX,  M(K*)  . 

XIII.  Calculate  DEXF  array  and  store  in  Table  5  (loop  6)  KK  =  2  to  HSTAT  +  1. 

A.  Set  DEXF  to  zero. 

B.  Accumulate  in  DEXF  sum  of  F2(LL,KK)  for  LL  =  1  to  KK  -  1  (loop  7). 

C.  Store  DEXF  in  Table  5  at  slot  KK  -  1. 

XIV.  Conditionally  recalculate  F2SUM  and  RDEX  in  NBW  section  storing  results 
in  Tables  3  and  4  (loops  23  and  25). 
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A.  List  heading  for  this  backwards  section. 

B.  (Loop  23)  Recalculate  F2SUM  and  RDEX  from  equations  (Dl)  and  (D2) 
for  L*  =  NSTAT  down  to  RSTAT  -  5  hy  using  stored  F2  values.  This 
uses  loop  24.  List  results  as  in  XII. E.  Store  results  in  TABLES 
3  and  4. 

C.  (Loop  25)  Continue  recalculation  of  F2SUM  and  RDEX  for  total  of  EB¥ 
gaps,  hut  this  time  use  proper  augmentation  and  deletion,  as  in 
loop  3.  This  uses  loops  26  and  27,  lists  results,  as  in  XII. E.,  and 
stores  results  in  TABLES  3  and  4. 

D.  Search  last  NBW  entries  in  TABLE  1  to  find  the  lowest  gap  with 
negative  entry  (loop  28). 

E.  Overwrite  TABLE  1  with  TABLE  3,  and  TABLE  2  with  TABLE  4  for  this 
gap  and  higher  gaps  (loop  29). 

XV.  List  TABLE  5  entries,  six  per  row  (loop  8). 

XVI.  List  tables  II  and  III  alternately  in  sets  of  twenty  gaps  per  page  start¬ 
ing  from  the  high  end  ( loop  9 ) . 

A.  Page  1  -  Last  twenty  gaps  -  table  II  (loop  21). 

B.  Page  2  -  Last  twenty  gaps  -  table  III  (loop  22). 

C.  Page  3  -  Next  to  last  twenty  gaps  -  table  II  (loop  21). 

D.  Page  4  -  Next  to  last  twenty  gaps  -  table  III  (loop  22). 

E.  Page  5  -  Remaining  gaps  -  table  II  (loop  2l). 

F.  Page  6  -  Remaining  gaps  -  table  III  (loop  22). 

XVTI.  Skip  page.  Return  to  step  VTI. 
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DICTIONARY  OP  SYMBOLS 


BIAS-  BIAS  IN  EXPONENT  NECESSARY  TO  KEEP  INTEGRAND  IN  RANGE 

DELE-  DELTA  E  OR  IP(L>-IP(K)  NEEDED  IN  VQEX  INTEGRATION 

DEXF  DEEXCITATION  FREQUENCY.  SEE  EQUATION  D3 

El-  IONIZATION  POTENTIAL  FOR  L  STATE  IN  VQEX  CALCULATION 

EN  ENERGY  LEVEL  FOR  STATE  IN  RECIPROCAL  CENTIMETERS 

F2  FUNCTION  2  DEFINED  BY  BY  EQUATION  D5 

GSIP-  GROUND  STATE  IONIZATION  POTENTIAL 

H-  PLANCKS  CONSTANT  IN  ERG-SECONDS 

INF-  INFINITY  FOR  VQEX  INTEGRATION 

IP-  IONIZATION  POTENTIAL  IN  ELECTRON  VOLTS 

K-  K  IN  ERGS  PER  DEGREE  KELVIN 

KSTAR-K*  THE  INDEX  OF  THE  UPPER  GAP  LIMIT 

KTHET-K  TIMES  THETA  IN  ELECTRON  VOLTS 

LAB  LABEL  FOR  A  STATE 

LKD-  LIMIT  ON  L-K  DIFFERENCE  TO  SHORTEN  EXECUTION  TIME 
LSTAR-L*  THE  INDEX  OF  THE  LOWER  GAP  LIMIT 
ME-  MASS  OF  ELECTRON  IN  GRAMS 

NBW-  NUMBER  OF  BACKWARD  STEPS  TO  TAKE  IN  GAP  CALCULATIONS 

NE  ELECTRON  DENSITY 

NK-  POPULATION  OF  STATE 

NSTAT-NUMBER  OF  STATES 

OMEGA-TOTAL  DEGENERACY  FOR  STATE 

PARI-  PARAMETER  1  DEFINED  IN  EQUATION  D6 

PAR2--PARAMETER  2  DEFINED  IN  EQUATION  D7 

PAR3-  PARAMETER  3  DEFINED  BY  EQUATION  D8 

SCAL-  SCALE  FACTOR  TO  COUNTERACT  BIAS  AFTER  INTEGRATION 

SIGZ-  SIGMA  ZERO  IN  CM  SQR-EV  SQR 

TAB-  TABLE  FOR  SUMMARY  PURPOSES 

TABLE  1  F2SUM  FROM  FORWARD  SWEEP.  SEE  EQUATIONS  D2  AND  D5 

TABLE  2  RDEX  FROM  FORWARD  SWEEP.  SEE  EQUATIONS  D1  AND  D4 

TABLE  3  F2SUM  RECALCULATED  IN  NBW  SECTION 

TABLE  4  RDEX  RECALCULATED  IN  NBW  SECTION 

TABLE  5  DEXF  DEFINED  IN  EQUATION  D3 
THETA-TEMPERATURE  IN  DEGREES  KELVIN 
WE  OMEGA  ELECTRON 


PROGRAM  LISTING 


SIBFTC  FAIN  DECK 

C  PROGRAM  FOR  STUDYING  RECOMBINATION  OF  AN  ATOMIC  ION 

DIMENSION  LABI  100)  ,£N.(  100)  ,  OMEGA  I  100  )  ,  F  2  (  99, 99  )  ,  NK  (  99  )  ,  I  P  (  99  ) 
COMMON  El, KTHET, BIAS, DELE 
DIMENSION  TAB(5,99) 

REAL  NE,NK,  IP,  INF,KTHET,K,ME,  KT,NKS 

PI=i. 141B92T 

K=1.38E-16 

H=6.625E-27 

SIG2=6. 53E-14 

ME=s9.1E-2B 

GS1P=31A07. 

NSTaT:‘47 
WE  =  2  * 

LKD=99 

READ  (5,20)  (LA8(J),EN(J),J=1,NSTAT) 

20  FORMAT  (6(15  ,F7.0)) 

READ  (5,30)  (OMEGA! J) ,J=1,NSTAT) 

30  FORMAT  (2AF3.0) 

WRITE  (6,40):  (  OMEGA  (  J  )  ,  J  =  1 ,  NSTAT  )  ,  GS  I  P  ,  WE 

WRITE  (6,40)  ( EN(J) ,J=1, NSTAT) 

40  FORMAT  (3X,10F9.1) 

N=NSTAT 

NP1=N+1 

NP2=N+2 

OMEGA! NPl )=10. 

LAB(NP1)=999 
EN(NP1)=GSIP 
ENINP2)=GSIP 
DO  1  J=1,NP2 

1  IPI  J).=  (GSlP-EN(J))/8067. 

WRITE  (6,50)  (  IP! J),J  =  1,NP2) 

50  FORMAT  (4H  IP  1P9E12.5) 

11  READ  (5,10)  NE, THETA, N6W 

10  FORMAT  (F10.0,F5. 0,415) 

NE=*NE»1.E  +  13 

WRITE  (6,60)  THETA, NE,NBW 

60  FORMAT  ( 1H1/5X , 7HTHETA=  ,  F 1 1 . 4 , 20X , 4HNE=  IPEIO.1,18) 

KT=‘K»THETA 
KTHET=‘KT/1.6E-12 

PAR1=1, 30035E-24*NE/KT*SIGZ/SQRT(KT)*PI/SQRT(ME) 
X=lta/KT)*(H/ME)/(2.*PI) 

PAR2=NE«NE*X*SQRT( X)/WE 

INF=.01*THETA 

DO  2  J^1,NP1 

2  NKl  J)=iPAR2*0MEGA(  J  )  »EXP  (  I P  (  J) /KTHET  ) 

WRITE  (6,70)  (NK( J) , J=1,NP1) 

70  FORMAT  ( 3H  XN  1P9E12.4) 

C  SECTION  TO  KEEP  INTEGRAND  IN  RANGE  FOR  LOWEST  GAP 
BIAS-0. 

EXPON*-IP(  D/KTHET 

12  IF(£XP(EXPON) )  13,13,14 

13  BIAS=BIAS+10. 

EXR0N=s£XP0N-10. 

GO  TO  12 

14  SCAL  =  l.E-iO*!EXP,(-BIAS) 

C  MAIN  LOUP  FOR  COMPUTING  RDEX  SUMS  FROM  LOWEST  TO  HIGHEST  GAP 
WRITE  (6,80) 
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80  FORMAT  (A3H1  GAP  L*  K*  F2SUM  RDEX  NK<K*) 

F2StjM:iO. 

RDEX=0. 

DO  3  LSTAR=1,N 

KSTAR»LSTARirl 

KMAX=KSTAR+LKD 

IF  IKMAX.GT.NPl)  KMAX=NP1 

LMirv  =  LSTAR~LKO 

IF  UMIN.LT*!)  LMIN=1 

C  LOOP  FOR  ADDING  AND  COMPUTING  LSTAR-K  TERMS 
El=*IPItSTAR> 

DO  4  KK=:KSTAR»KMAX 
DELt=IP(LSTAR)-IP(KK) 

PAR3=PAR1*$GAL*0MEGA(LSTAR) /OMEGA (KK)*EXP(OELE/KTHET) 

CALL  VQEXHDELE,  1NF,ANS) 

TERM=ANS»PAR3/DEL£**2 

D£L£=IP(LSTAR)-IP(KK+1) 

IF  1KK.GT.N3TAT)  GO  TO  15 
CALL  VQEX1{DELE, INFjANS) 

TERM=PAR3»ANS/0£LE**2“TERM 

15  TERM=AeS(TERM) 

F2{‘LSTAR,KK)=TERM 
F23UM^F2SUM-tTERM 
RDEX-RDEX+TERM*NKI KK) 

4  CONTINUE 

C  LOOP  FOR  SUBTRACTING  PROPER  L-KSTAR  TERMS 
LMI=LSTAR-1 

IF  iLMULT.l)  GO  TO  16 
DO  5  LL=LMIN,LM1 
T£RM=  F2(LL,LSTAR) 

F2SLM=tF2SUM“T£RM 
RDEX=RDEX-TERM*NK( LSTAR) 

5  CONTINUE 

16  TAB! 1,LSTARI=F2SUM 
TAB! 2|LSTAR1=RDEX 
NKS«=NKf  KSTAR) 

3  WRITE  (6t90>LA8(LSTAR) ,LAB(KSTAR) , L STAR , K STAR, F2 SUM, RDEX , NKS 

90  FORMAT  (219,^215,  1P4E22.4) 

C  SECTION  FOR  PUTTING  GUT  SUMMARY  TABLES 
DO  6  KK=2,NP1 
DEXF^O. 

DO  7  LL=2^KK 

7  DEXF=DEXF+F2( LL“1,KK) 

6  TAB!  5,KK-1)=JDEXF 

IF  INeW.:L£*0)  GO  TO  19 

C  NB.W  SECTION  FOR  RECALCULATION  OF  RDEX  AND  F2SUM  FROM  HIGH  GAP  END 

WRITE  J6,35)t  NBW 

35  FORMAT  {20H  BACKWARDS  SECTION  15, 6H  STEPS  ) 

DO  23  J=1^6 

I=NP1--J 

F2SUM^0. 

RD£X=0. 

JP=*I  +  i 

DC  24  KK  =  JP,INP1 
NKS=NK(KK) 

DO  24  L=1,I 
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24 


23 


26 


27 

25 

C 


28 


29 

19 

8 

31 


32 

33 


34 

21 

17 


TERM=F2(L,KK) 

F2Suy^F2SUM-tTERM 
ROEX=RDEX  +  TERM*IMKS 
TABl 3, 1 )=F2SUM 
NKS=KK( JP) 

TA614, I )=KDEX 

WRITE  16|90>  L A8 ( I ) » LAB ( JP ) » I » JP» F2SUM» RDEX » NKS 

DU  25  J=7»NBW 

I*NP1-J 

jr=i+2 


jp^  I+.1 
NKS=NK(JP) 

DD  26  L=lfl 
TERM=F2(Lf JP) 
F2SUM=»F2SUM  +  T£RM 
RDEX=RDEX+TERM*NKS 
DO  27  KK  =  JT,iNPl 
TERM=F2( JP.KK) 
F2SUM=tF2SUM-TERM 
RDEX=RDEX-TERM*NK(KK) 


TA6I3,  I)=F23UM 
TABl 4, I )=RDEX 

WRITE  (6,901  LABI  I) ,LAB(JP) ,1, JP,F2SUM,RDEX,NKS  ^ 

REPLACE  TABLE  1  AND  2  VALUES  WITH  BETTER  ONES  WHEN  AVAILABLE 


JT=‘G 

DO  28  J=1*NBW 
L=NP1-J 

IF  1TAB(  1,L) .GE.O. )  GO  TO  28 
JT=*L 

CONI INUE 

IF  IJT.LE.Ol  GO  TO  19 
DO  29  L=JT,N 


TABl 1,L)=TAB(3,L) 
TAB12,L)=TAB(4,L) 


DO  8  JA=UN3TAT,6 


JB^-JAtS 

WRITE  (6,311  JA, ( TAB( 5, J ) , J=JA, JB ) 
FORMAT  ( 3X, 15, 1P6E14. 5) 


KA=»NSTAT-19 

DO  9  JA= 1 ,N3TAT, 20 


KW=^2 

KB=*KA  +  20 

WRITE  (6,321  K6,KW 

FORMAT  (6H1  K=  I2,I5//2H  ) 

WHITE  (6,331  LA8(K8) ,NK{KB) 

FORMAT  ( 5X, 15, 28X, 1PE14. 2) 

DO  21  JB=1,20 
KB^KB-1 

IF  IKB.LT.ll  GO  TO  17 

AVIR  =  .5*(  IP(!KB1+IP(KB  +  1 )  1 

WRITE  (6,341  TAB( 1,KB) ,TAB(2,KB) , AVIP 

FORMAT  (■10X,ilP2E14.2,14X,OPF9.3) 

WRITE  16,331  LA8(KB) ,NK(KB) 

CONTINUE 


KW=3 

K8=*KA  +  20 
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VnRIJE  (6,32>  KBrKW 

WRilE  (6,33)  LA8(KB) ,NK(KB) 

DO  22  JB=i,20 
Ka=iKB-l 

IF  (KB^LE.O):  GO  TO  18 
WRITE  ( 6,34>  TABi 5,KB) ,TAB(2,KB) 
WRITE  (6,33)  LABIKB)  ,NK(KB)- 
22  CONTINUE 

9  KA=*KA-20 

18  KW^2 

WRITE  (6,32)t  Ke,KW 

GO  TC  11 

END 


SIBFTC  VGEXl 

SUBROUTINE  VQEXI  IXL,XU,ANS) 

X=^XL 

Y=XU 

0=1 Y-X)/512. 

XA=*X  +  D 

E=l.E-5*AbS(iF(  X)  ) 

12  IF  ( ABSIFIXA) ).LT.E)  GO  TO  11 

XA=«XA+D 

IF  (XA.LT^Y)  GO  TO  12 
II  Y=XA-D 

D=t V-X)/32. 

D2^L)  +  D 
N=  16 
T=X+D2 
S=F1T) 

DO  1  J=3,N 
T  =  T*tD2 

1  S=SfF(T) 

A=CF( X)+FtY>+2.*S)*D/3. 

T=X  +  D 
S^FI T) 

DO  2  J:i2,N 
T=T't-C2 

2  S=S+F(T) 

B=4. «S#D/3* 

R=AfB 

E=.G0001*ABS(‘R) 

IF  IE)  5,5,7 
7  DO  3  J=il,2 

C=.5«{ A+1.5<B) 

D2=iD 
D=*  5«0 
N^2tN 
T=X  +  D 
S=FI T) 

DO  4  K=2,N 
T=T+C2 

A  S^S^-FIT) 

B=4«*S*D/3« 


mm  33  m  fvj  ^ 


A*R-C 

W^AtiSie-C) 

IF  IVv-E)  5f5,3 

CONTINUE 

ANS^R 

RETURN 

END 


F 

FUNCTION  F(ARG)- 

COM^ICN  El  ,KTHET,3IAS,DELE 

REAL  KTHET 

E2=*ARG 

IF  IC6,GT,0^)  GO  TO  3 
C6:42*/3, 

EIfc2=El/E2 

DEE2=DELE/E2 

E2ES=E2/(EUE2) 

IF  nDELE  +  El).GT.E2)  GO  TO  2 

ECUATIGN  lA 

G*C6#E1E2+DEE2*( 1*-E1E2 )-DEE2»*2 
GO  TO  4 
EGUATICN  IB 

G=C6*( E1E2+DEE2*(1.-E1E2 )-DEE2**2) *SQRT ( ( 1 .+DELE/E1 ) *  1 1.-DEE2) ) 

F=E2Ei>*SQRTf  E2ES  )  *G*EXP  (  B I A  S-E2/KTHET  )  *E2*  1 .  E 10 

RETURN 

ENG 
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SAMPLE  DATA 


$cata 

600 

C 

601 

11178 

601 

11732 

701 

21766 

701 

21947 

602 

22589 

SOI 

25709 

801 

25792 

702 

26048 

504 

27010 

901 

27637 

901 

27681 

604 

28351 

1001 

28727 

1001 

28754 

1101 

29404 

1101 

29421 

1002 

29471 

903 

30043 

1301 

30170 

1202 

30198 

1103 

30495 

1501 

30566 

1402 

30581 

2  2 

4  4 

6  2 

2  4 

4  6 

2  14 

40  2 

4  10 

2  14 

2  4 

10  16 

6  10 

100  ICCO 

43 

100  5GOO  00 


502 

14499 

502 

14597 

700 

18535 

602 

22631 

800 

24317 

403 

24472 

702 

26069 

900 

26911 

.503 

26971 

802 

27816 

1000 

28300 

603 

28330 

902 

28832 

1100 

29130 

703 

29148 

1200 

29675 

1201 

29859 

1102 

29898 

1003 

30302 

1401 

30393 

1302 

30417 

1203 

30641 

1601 

30694 

2  4 

4  6 

2  14 

18  2 

4  10 

2  14 

14  6 

10  14 

6  10 

14  6 

10  14 

6 
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TABLE  I.  -  QUANTUM  NUMBERS,  IONIZATION  POTENTIALS,  AND  STATISTICAL  WEIGHTS  OF  ELECTRONIC  LEVELS  USED  IN  CALCULATIONS 


Electronic 

energy- 

level, 

K  or  L 

Principal 
and  secondary 
(orbital) 
quantum 
number, 
n,  1 

Spin  orbit 
quantum 
number, 

J 

State  binding 
energy  (ioni¬ 
zation  poten¬ 
tial), 

ev 

Total  statisti¬ 
cal  weight  of 
quantum  state 
(degeneracy), 

CDj^  =  Z(2J  +  1) 

Electronic 

energy 

level, 

K  or  L 

Principal 
and  secondary 
(orbital) 
quantum 
number, 
n,  1 

Spin  orbit 
quantum 
number, 

!  J 

State  binding 
energy  (ioni¬ 
zation  poten¬ 
tial), 

ev 

Total  statisti¬ 
cal  weight  of 
quantum  state 
(degeneracy), 

=  Z(2J  +  1) 

16p,p’ 

^16,1 

i  li 

2’ ■^2 

8.03  E-02 

6 

10s 

^10,0 

1 

2 

3.85  E-01 

2 

ief,f ' 

^12,3 

9.49  E-02 

14 

8d,d' 

1  ^8,2 

4.45  E-01 

10 

14d,d' 

®14,2 

2»  2 

1.02  E-01 

10 

9p’ 

9,1 

4.61  E-01 

4 

15p,p’ 

^15,1 

-,1- 
2  2 

1.04  E-01 

6 

9p 

9,1 

1 

2 

4.67  E-01 

2 

llf,f  ’ 

^11,3 

2-^,5^ 

2  2 

1.13  E-01 

14 

5g,g’ 

^5,4 

3i,4i 

2  2 

5.45  E-01 

18 

13d, d’ 

^13,2 

li,2i 

2  2 

1.22  E-01 

10 

5f,r' 

^5,3 

2i,3i 

2  2 

5.49  E-01 

14 

14P.P' 

^14,1 

2  2 

1.25  E-01 

6 

9s 

^9,0 

1 

2 

5.57  E-01 

2 

10f,f » 

^10,3 

pi  -zl 

22,  3g 

1.36  E-01 

14 

7d' 

7,2 

6,61  E-01 

6 

12d,d' 

^12,2 

ll,2l 

2  2 

1,49  E-01 

10 

7d 

7,2 

2 

6.64  E-01 

4 

13p,p’ 

^13,1 

2"  2 

1.53  E-01 

6 

8p’ 

8,1 

6.96  E-01 

4 

9f,f' 

®9,3 

2i,3l 

2  2 

1.69  E-01 

14 

8p 

8,1 

1 

2 

7.06  E-01 

2 

lld,d’ 

^11,2 

1~  2^ 

1.87  E-01 

10 

4f,f ' 

^4,3 

2i  si 

8.59  E-01 

14 

12p,p' 

^12,1 

2"  2 

1.92  E-01 

6 

0s 

^8,0 

1 

2 

8.78  E-01 

2 

0f,f ',123 

^8, 3, ^12,0 

2I  si  1 

2'  2 

2.14  E-01 

16 

6d' 

6,2 

1.07  E-00 

6 

10d,d' 

^10,  2 

li,2i 

2  2 

2.40  E-01 

10 

6d 

6,2 

li 

2 

1.09  E-00 

4 

lip'  . 

11,1 

2.46  E-01 

4 

7p' 

7,1 

1.17  E-00 

4 

lip 

11,1 

1 

2 

2.48  E-01 

2 

Vp 

7,1 

1 

2 

1.19  E-00 

2 

7f,f ' 

^7,3 

2i,3i 

2  2 

2.00  E-01 

14 

7s 

^7,0 

1 

2 

1.59  E-00 

2 

llB 

^11,0 

1 

2 

2.82  E-01 

2 

5d' 

5,2 

2.00  E-00 

6 

9d,d' 

^9,2 

li,2i 

2  2 

3.19  E-01 

10 

5d 

5,2 

2.09  E-00 

4 

lOp* 

10,1 

li 

2 

3.20  E-01 

4 

6p’ 

6,1 

jA 

2 

2.43  E-00 

4 

lOp 

10,1 

1 

2 

3,32  E-01 

2 

6p 

6,1 

1 

2 

2,50  E-00 

2 

6g,g'»h,h' 

a6,4,^6,5 

si  4i  4i  si 
*2,  =2 

3,79  E-01 

40 

6S 

^6,0 

1 

3.89  E-00 

2 

6f,f ' 

^6,3 

2  2 

3.81  E-01 

14 

Levels  degenerate  In  J  value. 


TABLE  III.  -  SUPERELASTIC  COLLISION  FREQUENCIES  FROM  ELECTRONIC  STATES  AND  DEEXCITATION  RATES 

ACROSS  ELECTRONIC  QUANTUM  GAPS 

[Free-electron  numUer  density,  10^^  cm”^;  free-electron  temperature,  5000^^  K.] 


Electronic 
state 
(47  to  1 
Inclusive), 
K 

Deexcitation 
■frequency 
from  state  K 
(superelastic), 
L=K-1 

E  "sek*’ 

L=1 

sec‘^ 

Deexcitation 
rate  across 
gap 

K*  to  L^^, 
^dex' 

cm‘^  sec"^ 

Popula¬ 
tion  of 
state  K, 

%’ 

-3 

cm 

Electronic 
state 
(47  to  1 
inclusive) 
K 

Deexcitation 
frequency 
from  state  K 
(superelastic) , 
L=K-1 

Yj  ^sek*' 

L=1 

sec“^ 

Deexcitation 
rate  across 
gap 

K*"  to  L"" , 
%ex> 

-3  -1 

cm  sec 

Popula¬ 
tion  of 
state  K, 

cm“^ 

K  47  16  ] 

4.32E  09 

K=24  6,3 

1.99E  10 

1  ?  3 

1.41E  13 

9.26E  22 

1.02E  10 

10,0 

3.61E  12 

1.14E  23 

2.87E  09 

14  2 

4.28E  12 

7.34E  22 

7.45E  09 

8,2 

1,51E  11 

4.71E  21 

1.65E  10 

15  1 

4.8GE  13 

3.98E  25 

4.49E  09 

9,1 

1.14E  12 

2.12E  22 

6.85E  09 

11  5 

3.42E  12 

6.07E  22 

1.07E  10 

9,1 

9.23E  12 

7.08E  22 

3.47E  09 

2.82E  12 

4.97E  22 

2,05E  11 

5.85E  21 

13  2 

7.eOE  09 

5,4 

2.91E  10 

2.18E  13 

1.94E  23 

2.60E  13 

7.61E  23 

14.1 

4.72E  09 

5,3 

2.94E  10 

2.76E  12 

4.31E  22 

9.89E  11 

4.75E  22 

10  3 

i.l3E  10 

9,0 

4.2eE  09 

1.79E  12 

3.37E  22 

5.88E  10 

2.13E  21 

12  2 

8.31E  09 

7,2 

1.63E  10 

1,78E  13 

1.64E  23 

6.17E  13 

I.OIE  24 

13, 1 

5.03E  09 

7,2 

l.lOE  10 

1.42E  12 

2.71E  22 

1.20E  11 

I.OIE  22 

9, 3 

1.22E  10 

8,1 

1.18E  10 

1.07E  12 

2.16E  22 

2.87E  12 

3.69E  22 

11,  2 

9.06E  09 

8,1 

6.04E  09 

1.05E  13 

1.05E  23 

3.53E  10 

1,67E  21 

12,1 

5.50E  09 

4,3 

6,03E  10 

e.57E  11 

1.81E  22 

2.64E  11 

1.63E  22 

12,0;8,3 

1.55E  10 

8,0 

9.01E  09 

5.54E  11 

1.41E  22 

2.00E  10 

I.OIE  21 

10, 2 

1.02E  10 

6,2 

4.39E  10 

6.16E  12 

7.01E  22 

1.66E  13 

7,29E  23 

11, 1 

4.16E  09 

6,2 

2.96E  10 

4,72E  13 

2.19E  23 

1.77E  10 

4,39E  21 

11,1 

2.09E  09 

7,1 

3.56E  10 

5.66E  11 

1.26E  22 

6.50E  11 

2.43E  22 

7,3 

1.57E  10 

7,1 

1.88E  10 

1.38E  13 

2.22E  23 

1.25E  09 

2.02E  20 

11, 0 

2.26E  09 

7,0 

4.75E  10 

3.80E  11 

l.lOE  22 

1.55E  10 

8,52E  20 

9.2 

1.23E  10 

5,2 

4.42E  11 

3.27E  12 

4.68E  22 

3,03E  12 

1.34E  24 

10, 1 

5.04E  09 

5,2 

3.03E  11 

2.20E  13 

1.31E  23 

7.10E  08 

2.43E  21 

10, 1 

2.54E  09 

6,1 

6.71E  11 

9.64E  11 

1.98E  22 

6.95E  10 

4.73E  22 

6,  4;  6,5 

5.65E  10 

6,1 

3.94E  11 

2,80E  13 

1.59E  24 

7.64E  07 

‘^1.32E  20 

6,0 

9.78E  12 

for  {6P-6B)  gap;  R^.-  -  I  (1. 32x10^0  cm^3  ^  6.6x10^^  cm"^  sec"^. 

min 
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aeronautical  and  space  activities  of  the  United  States  shall  he 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.” 

— National  Aeronautics  and  Space  Act  of  1958 
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